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1. INTRODUCTION 

Development of renewable energy sources (RES) offers a remarkable benefit to sustainable power 
generation in all electric power system with low environmental pollution emissions compared with 
conventional energy sources based power generation. The most important aspect of renewable power 
generation, well known as distributed generation (DG), is almost zero fuel cost. However, DG connected 
existing power system is raising concerns on complex problems of system operation and control. Basically, 
the power flow in transmission lines will be changed in direction significantly causing undesired conditions 
of power losses in the power system. Although DG can additionally increase energy efficiency and enhance 
the capability of power system, it can adversely impact on power losses due to the change of unidirectional 
power flow to bidirectional power flow without appropriate power management [1]. Typically, DG is 
operated at rated capacity to achieve economical investment aspect, whereas it may lead to high level of 
power losses and critical voltage profiles in the power system. Therefore, generation scheduling method 
plays a significant factor for power management eliminating power losses and complexity of system 
operation problems. 

Since power losses have been studied and presented in numerous aspects, some novel optimal power 
flows have been traditionally used in power losses reduction aspect [2]. Network reconfiguration was used 
for physically altering a network topology that aimed at change of the network power flow to reduce power 
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losses in a power system [3]-[4]. Transmission expansion planning aspect could be employed to improve 
network losses considering reactive power [5]. Interconnection planning of microgrid has been discussed for 
enhancing power losses and reliability [6]. In distribution system, power losses are a significant issue due to 
high R/X ratio. DG and capacitor locations were proposed for power losses reduction and ameliorating 
voltage profiles in distribution system [7]. In addition, switching devices allocation was analyzed in order to 
gain power loss reduction as well as service restoration [8]. Moreover, various methods were proposed in an 
area of DG units siting and sizing [9]-[11]. 

Although power loss minimization has been studied extensively, this aspect is still attracted the 
numerous research interests. Especially, the generation scheduling or power management based on power 
losses minimization has also been focused. For power management in distribution system, operational 
scheduling has been presented for overall benefit maximization and loss minimization based on electrical 
vehicles integrated system [12], [13]. Real-time energy management strategy was proposed considering 
operational constraints and power flow direction in microgrid [14], [15]. Demand response and renewable 
energy generation were considered for minimizing power loss using heuristic algorithm [16], [17]. Energy 
storage system (ESS) integrated microgrid was investigated to improve the loss minimization based on 
intermittent DG [18], [19]. 

This paper concentrates on an optimal generation scheduling method of power system considering 
maximum renewable energy harvesting and power loss minimization. The proposed method determines the 
significant variable of RES based power generation employing DG and ESS. However, DG accommodation 
and sizing cannot be adjusted in practical situation due to the capability constraints and economic benefit of a 
producer. Especially, DGs are always located at a specific areas that is uncontrollable factor. Although 
maximum DG outputs are expected for the producers, it can increase the power losses in the power system. 
Accordingly, DG output powers are managed within a proper way of maximum renewable energy harvesting. 
Then sufficient power will be supplied the loads while excessive power will be stored in ESS. The proposed 
method is investigated on IEEE standard 14-bus and 30-bus test systems via DIgSILENT Powerfactory 
software. Moreover, the proposed method is compared with an existing scheduling method to evaluate an 
effectiveness and confirm a superior performance. 


2. MATHEMATICAL MODEL 
The optimal generation scheduling problem for maximizing renewable energy harvesting from DG 
dispatch and minimizing power losses can be formulated and presented in (1) and (2) respectively. 


Maximize Í =max (P, DG dispatch ) 


(1) 


Minimize Js = min (Pross ine) a 


where f; and fz are the proposed objective functions. Ppg dispatch 18 the renewable energy harvesting from DG 
dispatch in the power system (MW). P, is the power loss in a transmission line (MW). 


oss line 


2.1. Renewable energy harvesting model 

In practical aspect, DGs are always operated at the maximum rated power output (Pp¢ output). 
This aspect may lead undesired conditions of power losses in the power system. On the other word, the 
utilities cannot directly control the power injected by DG. Therefore, this paper focuses on the renewable 
energy harvesting consisting of a DG dispatch (Ppa aisparcn) and an excessive power (Peiorage) between DG 
dispatch and maximum rated power output. The excessive power will be stored in ESS. The renewable 
energy harvesting function is given as follows: 


Pag dispatch P, DGoutput z; Pige (3) 


As the excessive power, which highly links with the power loss in ESS, is stored in ESS. 
The electric energy storage system (ESS) loss can therefore be detailed consisting of battery loss and 
converter loss [20]. Thus the ESS loss can be calculated as follows: 


P 


Loss ESS 


P 


Loss battery 


T Pris converter (4) 
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P =I} n x R patery (5) 


Loss battery battery 


P 


Loss converter T P, + (k% X Teog ) (6) 
where Pross battery and Pross converter are the battery loss and converter loss, respectively. Rpartery is the internal 
resistance of battery. Ipaiery is the charging current related to the amount of power stored in ESS (P sorage). Psp 
is the constant standby loss according to the power consumed by control platform, gate drivers, display, 
transducers and cooling fans. k% is the percentage of semiconductors and filter losses. 

However, this paper considers the direct relationship of the maximum renewable energy harvesting. 
Hence, ESS loss is particularly assumed to be a variable of Pyorage instead Of Ipanery. AS the Pstorage has been 
represented in (3) which significantly affects on the power loss of ESS. Consequently, the ESS loss can be 
assumed into a stored power in ESS and efficiency of ESS (7) that is provided as follows: 


P rige = P, DGoutput T P, DG dispatch (7) 


Possess = (1 E 7) P ige (8) 


2.2. Power loss in line model 

In order to obtain the power loss in line of the power system, the generalized power flow equation is 
employed in this paper to calculate the power losses. The power flow equation, deals with steady-state 
analysis related to real power and reactive power [21], can be expressed as follows: 


S,=F, + jQ, (9) 

P=d VV, Y, |cos(@, -O,-@,), i=1,2,....n (10) 
k=l 

0,=>'\V,||V, ||¥,.|sin(@, -O,-@,), i=1,2,....0 (11) 
k=l 


where S;, P; and Q; are the net apparent power, active power and reactive power injections to bus i, 
respectively. n is the total number of buses in the system. V; and V; are the voltage magnitudes at buses i and 
k, respectively. 0; and 6, are the voltage angles at buses i and k, respectively. Y is the admittance magnitude 
between buses i and k. aj, is the phase angle of admittance between buses i and k. This paper merely 
considers the power losses in lines with respect to active part from a branch conductance (g;,) between buses 
i and k that can be described as follows: 


P 


Loss line, ~ Sik [v +V — 2V;V, cos (8, — 4, | (12) 
2.3. Objective function formulation 
As Section 2.1 and 2.2, the maximum renewable energy harvesting is to maximize power dispatch 
of DG or minimize excessive power which can be represent minimum the power loss in ESS. The 
combination of Equations (1) and (2) is the objective function of the proposed method. Therefore, it can be 
expressed as follows: 
NI Nst 
Min Pe DN line,i + >a ESS, j (13) 
j=l 


i=l 


where Pross iinei 18 the power loss in line i, Pross gssj is the ESS loss at location j. NI and Nst are the total 
number of lines and the total energy storage locations, respectively. 
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2.4. Operational constraints 
2.4.1. Power flow constraint 

The power flows across each line from any two buses i and j must be remained within the maximum 
limit by maximum current capacity of each line. The power flow constraint is given as follows: 


| ee es (14) 


where J;j is the current in the line between buses i and j. J, is the maximum current capacity of the line 


between buses i and j. 


2.4.2. Generator constraints 

The generators in the system must be operated under the rated active and reactive power limits 
associated with the bus voltage. The voltage level must be also within maximum and minimum limits. 
The generator constraints can be provided as follows: 


P< P, < PP” (15) 
Or < Qy SOR™ (16) 
Vee (17) 


where Py and Q; are the active and reactive power injection at generator bus N, respectively. P™ and Q™ 
are the maximum active and reactive powers of the generator N, respectively. P™ and Q"" are the minimum 
active and reactive powers of the generator N, respectively. Vy is the bus voltage where a generator 
connected at bus N. V,"" and V,"" are the maximum and minimum operating voltage limits of a generator 
bus, respectively. 


2.4.3. Renewable distributed generation constraint 
The renewable DG constraint is only considered the power output limit. The limit can be given as 
follows: 


0 < Poon S Poon (18) 
where Ppg,v is the active power dispatch from DG to bus N. P" is the maximum rated active power of each 
DG at bus N. 


2.4.4. Load constraints 

The general loads are distributed among the system while the loads should be operated in the 
voltage constraints as given in (19). Additionally, the load must be also operated under the voltage deviation 
(VD) limit. VD is represented as the difference of voltage between maximum and minimum voltage limits. 
The VD can be expressed as given in (19) and (20). 


Vo <V, SVP“ N =1,...,n. bus no. (19) 


VD, = V™™* —yVy™ i =1,...,m. scinarios no. 20) 


where Vy is the bus voltage where the load is connected at bus N. Vand Vare the maximum 
and minimum operating bus voltage limits, respectively. V™ and V™ are the maximum and minimum of 
system voltage at scenario i, respectively. 
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3. PROPOSED METHOD 

This proposed method applies PSO algorithm to explore the optimal solutions of generation 
scheduling for maximizing renewable energy harvesting and minimizing power losses in accordance with 
Section 2. The PSO algorithm is initialized with a set of random positions of particles (solutions) [22], [23]. 
Each particle will be executed to obtain the solution according to the objective function. Every iteration, 
the particles will be updated by following two important values. The first one is the best solution (Pj.s,) that 
links with cognitive factor. The other important value is tracked by the particle swarm optimizer in 
accordance with social factor. This value is a global best (gpest). The particles are updated according to the 
Poest aNd Zpest for obtaining the new positions and velocities using (21) and (22). 


k 


i 


kk ko ky ko 
Vi FV TON Presti —% Ch | E besti — X 


(21) 


i 


— 


k+l _ „k k+l 
X =X +v; (22) 


where v‘ is the velocity of a particle i at iteration k. x‘ is the position of a particle i at iteration k. Prag 1S 
the personal best position of the particles. gý „„; is the global best position of the particles. c; and cz are 


positive acceleration constants. r; and r, are the random values from uniform distribution. 

Although the PSO algorithm is initialized with the random positions of particles x to search the best 
position for obtaining the optimal solutions with respect to the proposed objective function, the swam can be 
determined as the multidimensional variables according to a number of designated variables in the power 
system. Therefore, the sets of particles (X) and velocities (V) are associated with a number of variables (n) 
and number of particles (m) that can be presented in (23) and (24). This paper considers 3 variables 
consisting of the conventional power generation (Py), voltage at generation bus (Vy) and renewable 
distributed generation (Pg aisparcn). The flowchart of the proposed method is represented in Figure 1. 


Initalize particles in swam and velocities 
under the constraints 


For iteration’ =1 


For particle i= 1 


Execute particle ito obtain the objective functions 
Evaate and set Daw 


Next particle #1 


Maximum particle 
isreached 


Update the positions 
and velocities 
according to 
Zor and pwe 


Evaluate gav: and Pweaccarding to 
1. Max. Renewabk energy harvesting 
2. Min. Power losses 


Maximum iteration 
isreached 


End o process 


Figure 1. Flowchart of PSO based proposed method 
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4. CASES FOR COMPARISON 

In this section, two cases are performed for comparative study in order to exploit the capability of 
the proposed method. According to the objective function in Section 2, the study was designed in 
DIgSILENT Program Language (DPL). Furthermore, all cases are operated under the designated operational 
constraints as described in Section 2. The performed cases can be detailed as follows: 


4.1. Line loss minimization (case1) 

The generation units (conventional generation units and DG units) in the system will be searched to 
receive the optimal power dispatches considering solely power losses minimization. The power loss 
minimizing is considered the power losses in transmission lines under the operational constraints. 
Additionally, the search condition is regarded to voltage deviation. Therefore, the optimal solutions of case 1 
(xı) can be defined according to the designated variables including conventional generation power (Py), 


voltage level (Vy), and DG power ( P% ) that can be provided as follows: 


DG dispatch 
e P V. p% 1 E (2 5 ) 
% S| in VN !pGdispatch 


pr! e [0, Pyg joa o 


DG dispatch 


4.2. Maximum DG output dispatch (case2) 

The maximum renewable energy harvesting is the main objective function of this case. The method 
will be explored the optimal point of conventional generations and maximizing renewable energy harvesting 
from DG units. The search space of DG power is consistent and recalled from case1. However, the optimal 
renewable DG dispatch must to search within the DG limit and the previous DG power dispatch. The optimal 
solution of case 2 (x2) can be provided as follows: 


T 
a opt.2 
Ag Es w Vy PpGdispatch ] (27) 
opt.2 opt.1 
Poo dispatch € LP, DG dispatch? P, DG output ] (28) 


5. SIMULATION RESULTS AND DISCUSSION 
5.1. Test systems description 

The proposed method is tested in comparative study with the performed cases under IEEE 
standard 14-bus and 30-bus test systems. The generation units consist of the conventional generation units 
and renewable DG units according to the test systems. The ESSs are installed in each location of renewable 
DG unit to collect the excessive power. The available DG power output at each location is comprised of 
power energy available of DG dispatch and stored excessive power in ESS. The components data of the 
14-bus and 30-bus test systems are detailed in Table 1 and Table 2, respectively. To assess the ESS losses, 
the efficiency of ESS is assumed to be 90% for all installed locations. 
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Table 1. Detail of IEEE Standard 14-bus Test System Table 2. Detail of IEEE Standard 30-bus Test System 


Cap. Cap. 
No. Type Bus [MW] No. Type Bus [MW] 
1 Conventional gen. unit 1 1 750 1 Conventional gen. unit 1 1 200 
2 Conventional gen. unit 2 2 600 2 Conventional gen. unit 2 2 150 
. : 3 Conventional gen. unit 3 5 150 
3 Conventional gen. unit 3 3 400 4 Conventional gen. unit 4 8 80 
4 Renewable DG unit 1 12 100 5 Conventional gen. unit 5 11 50 
5 Renewable DG unit 2 10 100 6 Conventional gen. unit 6 13 50 
6 Renewable DG unit 3 9 100 T Renewable DG unit 1 5 50 
S unit 1 2 8 Renewable DG unit 2 3 50 
1 ESS unit l 9 Renewable DG unit 3 9 50 
8 ESS unit 2 10 10 ESS unit 1 5 
9 ESS unit 3 9 11 ESS unit 2 3 
12 ESS unit 3 9 


5.2. Simulation Results 

This paper considered the variation of each generated DG power output in a particular day as shown 
in Figure 2. In general, DGs are always generated at maximum rated capacity based on renewable energy 
resources in time interval of a day. The power dispatching of DGs is respected to power loss in line 
and power loss in energy storing process. 


110 ç 110 p 
100 Ẹ 


Active Power (MW) 
Active Power (MW) 


1 5 9 13 17 21 25 29 33 37 41 45 49 
1 5 9 13 17 21 25 29 33 37 41 45 49 Pee 
eer Scenarios in a day 
Scenarios in a day 


(a) (b) 


110 g 


Active Power (MW) 


1 5 9 13 17 21 25 29 33 37 41 45 49 
Scenarios in a day 


(c) 


Figure 2. Power generation of renewable distributed generation in a particular day (a) DG unit 1 
(b) DG unit 2 (c) DG unit 3 


5.2.1. IEEE standard 14-bus system 

The IEEE standard 14-bus test system was determined by 48 operating conditions in accordance 
with half an hour forecasting of load demand during a particular day in order to investigate the capability of 
the proposed method. The optimal generation scheduling was obtained by the proposed method compared 
with the performed cases as illustrated in Figure 3. Figure 3 demonstrated the load demand (MW), 
conventional generation (MW), maximum DG power output (MW) and DG dispatch (MW) in a particular 
day. The load demands were supplied by conventional generation together with DG to fulfill the mismatch 
power. The stored excessive powers in ESS were represented between the area of DG output curve and 
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optimal DG dispatch. Since the excessive power influences the power loss in stored power process, thus the 
level of power loss in power storing process is presented by this area. The simulation results were 
summarized in Table 3. 


1000 1000 
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—+— DG Output —#— DG dispatch —+— DG output —a— DG dispatch 
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Active Power (MW) 
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Scenarios in a day 
(c) 


Figure 3. Optimal generation scheduling and load demand (a) case 1 (b) case 2 (c) Proposed method 


Table 3. The IEEE 14-bus System Simulation Results 


Parameters Casel Case2 Proposed method 
DG Harvesting (MW) 6,840.9 7,879.6 7,509.9 
Loss in lines (MW) 228.9 313.7 226.7 
ESS Loss (MW) 210.9 107 143.8 
Total Power loss (MW) 439.8 420.7 370.5 
Max. VD (p.u.) 0.0999 0.118 0.1 


In case 1, the simulation results were demonstrated at Figure 3(a). Some scenarios in a particular day 
were focused. The load demand at scenario 5 was 614.057 MW. The conventional generations and DGs were 
optimally generated at 478.678 MW and 139.219 MW, respectively. The total excessive powers were stored 
77.780 MW as illustrated in Figure 4(a). Hence, the total power losses were 11.008 MW as illustrated in 
Figure 4(b). To make the point clear, scenario 25 was also focused. The load demand was 744.880 MW, 
whereas the conventional generations and DGs were optimally generated at 567.691MW and 182.571 MW, 
respectively. The total excessive powers were stored 83.4294 MW in ESS resulting the total power losses 
13.116 MW as illustrated in Figure 4(a) and Figure 4(b), respectively. 

In case 2, the scenarios in a particular day were recalled from the previous case. The optimal 
conventional generations were totally generated at 450.806 MW, and DGs were optimally dispatched at 
172.661 MW as illustrated in Figure 3(b). The total excessive powers were stored 44.338 MW in the ESS. 
The total power losses were 13.239 MW in scenario 5 as illustrated in Figure 4(b). In scenario 25, DGs were 
generated at 540.515 MW and the conventional generations were generated at 211.262 MW to meet the 
optimal solutions as illustrated in Figure 3(b). The excessive powers were stored 54.738 MW leading the 
total power losses 11.693 MW as illustrated in Figure 4(a) and Figure 4(b). 
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In proposed method, maximum renewable energy harvesting and minimum power loss were to be 
the objective functions. The conventional generations and DGs were optimally generated at 446.723 MW and 
172.318 MW to supply 614.057 MW of load demand in scenario 5 as illustrated in Figure 3(c). The total 
excessive powers were stored 44.681 MW and the total power losses were 8.8106 MW as illustrated in 
Figure 4(a) and Figure 4(b). The conventional generation and DG were optimally generated at 529.986 MW 
and 210.093 MW in scenario 25 to supply the load demand 744.880 MW as illustrated in Figure 3(c). The 
total were stored 55.907 MW resulting the total power losses 8.703 MW as illustrated in Figure 4(a) and 
Figure 4(b), respectively.excessive powers. 
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Figure 4. (a) The stored excessive power in ESS (b) Total power losses at the optimal schedule condition 


The voltage profiles were extensively investigated by voltage deviation (VD). Figure 5 illustrated 
the level of VD depending on the proportion of DG in system for all cases. The worst case was obviously 
indicated in case 2 which the objective function of this case was maximum renewable energy harvesting. Due 
to the voltage constraints were addressed according to Section 2, hence the system voltages must be 
improved. Consequently, the proposed method was controlled the VD within 0.1 p.u. during maximum 
renewable energy harvesting was considered. This implies that all buses in the system were still maintained 
within the marginal constraints where the minimum voltage is not reached below 0.95 p.u. 

The rest of this simulation study aimed at the convergence rate of the proposed method. 
The convergence rate of proposed method is depicted in Figure 6(a). The power losses are minimized to 6 
MW for scenario 46 which was approximately identified at iteration 6. However, the proposed method was 
applied by changing the search space, which it did not consistent with case 1 and case 2, for PSO algorithm 
to confirm its capability. At 6 MW of total power losses, the method was approximately converged to the 
minimum power loss after 95 iterations as illustrated in Figure 6(b). 
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Figure 5. DG generations and voltage deviations (a) case | (b) case 2 (c) Proposed method 
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Figure 6. (a) The PSO convergence rate of case 3 for 14 bus simulation result, (b) The PSO convergence rate 
of scenario 46 with changing search space for 14-bus simulation result 


5.2.2. IEEE standard 30-bus system 

To confirm the capability of the proposed method, the IEEE standard 30-bus test system was 
employed in this section. The performed cases were still recalled for investigating the proposed method. The 
simulation results demonstrated the superior effectiveness and performance of the proposed method to 
achieve the optimal solution based on the objective functions in Section 2. The simulation results are 
summarized in Table 4. 

The total DG dispatch over the operating scenarios in a particular day of each case studies was 
optimally dispatched at 436.721 MW in case 1, 588.144 MW in case 2, and 573.365 MW in proposed 
method case. Increasing of DG influenced the decrease in conventional generations and reducing of ESS loss. 
In cases 1 and 2, the generations were nevertheless dispatched with undesired solution because the total 
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losses were still high as provided in Table 4, and VD of the buses were reached the voltage constraints in 
some operating scenarios. The ESS loss was occurred as 39.509 MW in case 1, 30.260 MW in case 2, and 
31.737 MW in case 3. Hence, the total power losses for optimal generation scheduling were obtained as 
81.595 MW in case 1, 75.565 MW in case 2, and 65.847 MW in proposed method case as provided 
in Table 4. 


Table 4. The IEEE 30-bus System Simulation Results 


Parameters Casel Case2 Proposed method 
DG Harvesting (MW) 436.721 588.144 573.365 
Loss in lines (MW) 42.086 45.305 43.061 
ESS Loss (MW) 39.509 30.260 31.737 
Total Power loss (MW) 81.595 75.565 74.798 
Max. VD (p.u.) 0.1 0.108 0.105 
Convergence (Iteration) 10 8 5 


5.3. Discussion 

The power delivered from the generation units to the consumer points is always accompanied with 
power losses. The variations of generations can directly result the power losses in the system. The trending of 
power losses does not only depend on the conventional generation, but available of renewable DG may also 
increase power losses in the system and raise the complexity of power management. Hence, power losses 
should be dealt with the combination of generation types in the power system. 

The simulation results of comparative study demonstrate the different combination of power losses 
at the same operating scenarios. Although the ESS losses is obviously reduced in case 2, however, this 
reduction must be still included power losses in lines because the maximum dispatch of DG can result the 
power losses in lines especially long-distant transmission lines. 

The proposed method case is to minimizing total power losses and maximizing renewable energy 
harvesting. Although the proposed method is allowed the higher the ESS losses than the ESS losses in case 2, 
the total power losses in lines are reduced by increasing some the ESS losses in a way of maximum 
renewable energy harvesting. In Table 3 and 4, the simulation results demonstrate that the total power losses 
of the proposed method is significantly lower than the cases for comparison. 

The reduction of power losses are related to the optimal power for storage and minimize power loss 
in line. Figure 7 reveals the benefits of renewable energy harvesting and power loss reduction. This is the 
main advantage of the proposed method including two sections: the reduction of power losses have been 
improved and increased of renewable energy harvesting. 


—a— Increasing penetration of DG 


—— Power losses reduction 


Avtive Power (MW) 


Scinarios in a day 


Figure 7. Relationship of power losses reduction and increasing of DG dispatch 


5. CONCLUSION 

In this paper, an optimal generation scheduling has been investigated in the power system. 
The proposed method was executed considering maximum renewable energy harvesting and minimizing 
power losses. The optimal solutions of the proposed method were identified and obtained using PSO 
algorithm. The two cases for comparison were performed to exploit the capability and efficiency of the 
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proposed method. The simulation results demonstrated the effectiveness and performance of the proposed 
method to achieve the optimal solutions for generation scheduling especially with maximum renewable 
energy harvesting and minimizing power losses. The power losses were evidently decreased according to the 
optimal storing power in ESS and minimizing line losses with maximum renewable energy harvesting. 
Additionally, the maximum renewable energy harvesting is significantly influenced the decrease in 
conventional generations and reducing ESS losses. 
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